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Abstract 



The spherical gravitational collapse of a compact packet consisting of per- 
fect fluid is studied. The spacetime outside the fluid packet is described by 
the out-going Vaidya radiation fluid. It is found that when the collapse has 
continuous self-similarity the formation of black holes always starts with zero 
mass, and when the collapse has no self-similarity, the formation of black holes 
always starts with a finite non-zero mass. The packet is usually accompanied 
by a thin matter shell. The effects of the shell on the collapse are also studied. 

04.20.Jb, 04.40.Dg, 97.60.Lf. 
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I. INTRODUCTION 



Critical phenomena in gravitational collapse have attracted much attention |l] since the 
pioneering work of Choptuik [0]. From the known results the following emerges 0: In 
general critical collapse of a fixed matter field can be divided into three different classes 
according to the self-similarities that the critical solution possesses. If the critical solution 
has no self-similarity, continuous or discrete, the formation of black holes always starts with 
a mass gap (Type I collapse), otherwise it will start with zero mass (Type II collapse), and 
the mass of black holes takes the scaling form Mbh oc (P — P*) 7 , where P characterizes 
the strength of the initial data. In the latter case, the collapse can be further divided into 
two subclasses according to whether the critical solution has continuous self-similarity (CSS) 
or discrete self-similarity (DSS). Because of this difference, the exponent 7 is usually also 
different. Whether the critical solution is CSS, DSS, or none of them, depending on both 
the matter field and the regions of the initial data space The co-existence of Type I 
and Type II collapse was first found in the SU(2) Einstein- Yang-Mills case ||, and later 
extended to both the Einstein-scalar case || and the Einstein-Skyme case ||, while the 
co-existence of CSS and DSS critical solutions was found in the Brans-Dicke theory ]7|. The 
uniqueness of the exponent in Type II collapse is well understood in terms of perturbations 
||, and is closely related to the fact that the critical solution has only one unstable mode. 
This property now is considered as the main criterion for a solution to be critical Q. 

While the uniqueness of the exponent 7 crucially depends on the numbers of the unstable 
modes of the critical solution, that whether or not the formation of black holes starts with 
a mass gap seemingly only depends on whether the spacetime has self-similarity or not. 
Thus, even the collapse is not critical, if a spacetime has CSS or DSS, the formation of black 
holes may still turn on with zero mass. If this speculation is correct, it may have profound 
physical implications. For example, if Nature forbids the formation of zero-mass black holes, 
which are essentially naked singularities ffl, it means that Nature forbids solutions with self- 



similarity P]. To study this problem in its generality term, it is found difficult. In flOfl, which 



2 



will be referred as Paper I, gravitational collapse of massless scalar field and radiation fluid 
is studied, and it was found that when solutions have CSS, the formation of black holes 
indeed starts with zero-mass, while when solutions have no self-similarity it starts with a 
mass gap. 

In this Letter, we shall generalize the studies given in Paper I to the case of perfect fluid 
with the equation of state p = kp, where p is the energy density of the fluid, p the pressure, 
and k an arbitrary constant, subjected to < k < 1. We shall show that the emerging 
results are consistent with the ones obtained in Paper I. Specifically, we shall present two 
classes of exact solutions to the Einstein field equations that represent spherical gravitational 
collapse of perfect fluid, one has CSS, and the other has neither CSS nor DSS. It is found 
that such formed black holes usually do not have finite masses. To remedy this shortage, 
we shall cut the spacetime along a time-like hypersurface, say, r = r (t), and then join the 
internal region r < ro(t) with an asymptotically flat out-going Vaidya radiation fluid, using 
Israel's method JTTJ. It turns out that in general such a junction is possible only when a 
thin matter shell is present on the joining hypersurface [TJ|]. Thus, the finally resulting 



spacetimes will represent the collapse of a compact packet of perfect fluid plus a thin matter 
shell. The effects of the thin shell on the collapse are also studied. It should be noted that 
by properly choosing the solution in the region r > ro(t), in principle one can make the 
thin shell disappear, although in this Letter we shall not consider such possibilities. The 
notations will closely follow the ones used in Paper I. 

II. EXACT SOLUTIONS REPRESENTING GRAVITATIONAL COLLAPSE OF 

PERFECT FLUID 

In this section, we shall present two classes of solutions to the Einstein field equations, 

R P v - -^QtiuR = (p + p)u^u u - pg^, 

where is the four-velocity of the perfect fluid considered. The general metric of spherically 
symmetric spacetimes that are conformally flat is given by |L0 



ds 2 = G(t, r) [dt 2 - h 2 (t, r) (dr 2 + r 2 dQ 2 )] , (1) 

where dVl 2 = d6 2 + sin 2 6*<iy9 2 , {x^} = {t, r, 6 1 , </?} (// = 0,1,2,3) are the usual spherical 
coordinates, and 

h(t,r) = \ (2) 

with 7^ 0. The Friedmann- Robertson- Walker (FRW) metric corresponds to G(t,r) = 
G(t) and f\{t) = Const. The corresponding Einstein field equations are given by Eqs.(2.20) 
- (2.23) in Paper I. Integrating those equations, we find two classes of solutions. In the 
following, we shall present them separately. 

a) Class A Solutions: The first class of the solutions is given by 

G{t,r) = (1 -Pt) 2€ , h{r) = l, 

p = kp = 3keP 2 (1 - Pt)- 2 ^ +1) , u„ = (1 - Ptyt 51, (3) 

where P is a constant and characterizes the strength of the solutions (See the discussions 
given below), and £ = 2/(1 + 3k). This class of solutions is actually the FRW solutions 
and has CSS symmetry ||13|| . However, in this Letter we shall study them in the context of 
gravitational collapse. 

To study the physical properties of these solutions, following Paper I we consider the 
following physical quantities, 

p tl p2 3 

m f(t,r) = ^(l + R a R p g^ 
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n = R aPlS R - _ pt ^ (1H) , (4) 

where R is the physical radius of the two sphere t,r = Const., and m^(t,r) is the local 
mass function |fH| . From Eq.(^|) we can see that the spacetime is singular on the space- 
like hypersurface t = P~ x . The nature of the singularity depends on the signature of the 
parameter P. When P < 0, it is naked, and the corresponding solutions represent white 



holes [pT0| . When P = 0, the singularity disappears and the corresponding spacetime is 



Minkowski. When P > 0, the singularity hides behind the apparent horizon, which locates 
on the hypersurface, 

_ l-Pt 

r = r AH = , (5) 

with tah being a solution of the equation R !a R j/3 g a/3 = 0. Thus, in the latter case the 
solutions represent the formation of black holes due to the gravitational collapse of the fluid. 
The corresponding Penrose diagram is similar to that given by Fig. 1(a) in Paper I. Note that, 
although the spacetime singularity is always space-like, the nature of the apparent horizon 
depends on the choice of the parameter k. In fact, when 1/3 < A; < 1, it is space-like; when 
k = 1/3, it is null; and when < k < 1/3, it is time-like. Substituting Eq.(|5]) into Eq.(|), 
we find that m AH (t,r AH ) = {P0^ r AH- Thus, as r AH — > +oo, we have m AH — > +oo. That 
is, the total mass of the black hole is infinitely large. To get a physically reasonable model, 
one way is to cut the spacetime along a time-like hypersurface, say, r = ro(t), and then join 
the part r < r (t) with one that is asymptotically flat |J. We shall consider such junctions 
in the next section. 

(3) Class B Solutions: The second class of solutions are given by 

G(t, r) = sinh 25 [2af -1 (t ~ et )} > h ( r ) = ~ ^ 2 y\ (6) 



where e = ±1, £ is defined as in Eq.(y), £ an d a{= \/—fi) are constants. Introducing a new 
radial coordinate f by df = h(r)dr, the corresponding metric can be written in the form 



ds 2 = sinh 2 «[2r 1 (to - et)} jrft 2 - dr 2 - sinh ^( 2r ) ^j 



(7) 



Note that in writing the above equation we had, without loss of generality, chosen a = 1, 
and dropped the bar from f. The energy density and four-velocity of the fluid are given, 
respectively, by 

p = kp = 12ksm~ 2 ^[2(t - et)/£], UfM = sm^[2{t - et)/^, (8) 
while the relevant physical quantities are given by 



m f (r, *) = 7 sinh 3 (2r) sinli^ 2 [2£ _1 (£ - et) 

ft = 288 (l + e 2 ) r 2 sinh" 4(f+1) \2C\to - et) 



(9) 



The apparent horizon now is located at 



fAH = i (t -et). 



(10) 



From Eq.(^|) we can see that the solutions are singular on the hypersurface t = et^. When 
e = — 1 it can be shown that the corresponding solutions represent cosmo logical models with 
a naked singularity at the initial time t = —to, while when e = +1 the singularity is hidden 



behind the apparent horizon given by Eq.(10), and the solutions represent the formation of 
black holes due to the collapse of the fluid. In the latter case the total mass of black holes is 
also infinitely large. To remedy this shortage, in the next section we shall make "surgery" 
to this spacetime, so that the resulting black holes have finite masses. 



III. MATCHING THE SOLUTIONS WITH OUTGOING VAIDYA SOLUTION 

In order to have the black hole mass finite, we shall first cut the spacetimes represented 
by the solutions given by Eqs.(Q), and (0) along a time-like hypersurface, and then join the 
internal part with the out-going Vaidya radiation fluid. In the present two cases since the 
perfect fluid is comoving, the hypersurface can be chosen as r = r = Const. Thus, the 
metric in the whole spacetime can be written in the form 

[ A(t, r) 2 dt 2 - B{t } r) 2 dr 2 - C(t, r) W, (r < r ), 
ds 2 = < (11) 

{ i 1 - 22 R a ) ^ + 2dvdR2 ~ R2dQ2 > ( r ^ r °)' 

where the functions A(t,r), B(t,r) and C(t,r) can be read off from Eqs. ([!]), (||) and (^). 
On the hypersurface r = r the metric reduces to 

ds 2 \ r=ro = g ab di a di b = dr 2 - R(r) 2 dn 2 , (12) 

where £ a = {r, 6, <£>} are the intrinsic coordinates of the surface, and r is defined by 
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dr 2 = A 2 (t,r )dt 2 



1 - J ^ 2 + 2dvdi2, 



(13) 



where v and i? are functions of r on the surface, and R(r) = C(t,r ), M(r) = m{v{r)). 
The extrinsic curvature oon the two sides of the surface defined by 



K. 



ab 



-nt 



d 2 x± ± a dx±dx± 



has the following non-vanishing components |y| 



(14) 



k: 



B ' 
i; cM(r) 



sin" 2 OK 



R2 , K+=sm~ 2 6K^ = R Iv 1 



2M(r) 



(15) 



where i = dt/dr, () )ft = d( )/dx^ and are the normal vectors defined in the two faces of 



the surface. Using the expression |TT 



[K ab \- - g ab [K) =-8 



7TTn. 



(16) 



we can calculate the surface energy-momentum tensor r ab , where [-/£<,&] ~ = K^ b —K~ b , [K)~ = 
g ab [K ab ]~ , and g ab can be read off from Eq. fll2"D . Inserting Eq.fllSl) into the above equation, 
we find that r ab can be written in the form 



Tab = (TW a W b + 7] (9 a 9 b + (f) a (j) b ) , 



(17) 



where w a , 9 a and a are unit vectors defined on the surface, given respectively by w a = 
51, 9 a — R5 e a1 (p a = R sin 95%, and a can be interpreted as the surface energy density, r\ the 



tangential pressure, provided that they satisfy certain energy conditions ||17|1 . In the present 
case a and r\ are given by 



rr -K:lR-- + J'(r ) 



A-nR 



v 



V 



±-{v 2 -2vR-2vf(r ) + l}, 



(18) 



16ttRv 

where J(r) = r for Class A solutions, and J(r) = sinh(2r)/2 for Class B solutions, and 
a prime denotes the ordinary differentiation with respect to the indicated argument. Note 



that in writing Eq.flT8|) we had used Eq. (|13|) , from which it is found that the total mass 
of the collapsing ball, which includes the contribution from both the fluid and the shell, is 
given by 

M(r) = ^{v 2 + 2vR-l). (19) 

To fix the spacetime outside the shell we need to give the equation of state of the shell. In 
order to minimize the effects of the shell on the collapse, in the following we shall consider 
the case rj — 0, which reads 

v 2 - 2vR - 2J'(r )v + 1 = 0. (20) 
To solve the above equation, let us consider the two classes of solutions separately. 



A. Class A Solutions 



In this case, it can be shown that Eq. (p0|) has the first integral 



v(r) = X " 2( "°" 1) ^ , (21) 
V ; x-2v Ro ' V ' 

where R(r) = Rqx^, Rq = roP^ 1 , x ee [(£ + l)(r — r)]«+i, and v and r are integration 



constants. Substituting the above expressions into Eq . (|T9|) , we find that 

Mix) = [x _^lf_ 1 1)Ro? {( 2 " + 2(e - mv - l)Rox + 4^ (1 - «b)J2g} . (22) 

At the moment r = t^h (or x = xah — £,Ro), the shell collapses inside the apparent horizon. 
Consequently, the total mass of the formed black hole is given by 

M BH = M(x AH ) = ^_ 2( °^_ 1)]2 {^(2 - + 2(f - 1)(2« - 1) + 4u (l - v )} , (23) 

which is finite and can be positive by properly choosing the parameter vq for any given £. 
The contribution of the fluid and the thin shell to the black hole mass is given, respectively, 



S 



byQ 

f — f / \ S r D £ 
m M = m AH\' I ~AH) = P > 

mg^ = ^WM^) = f^^^- (24) 

From the above equations we can see that all the three masses are proportional to P, the 
parameter that characterizes the strength of the initial data of the collapsing ball. Thus, 
when the initial data is very weak (P — ► 0), the mass of the formed black hole is very 
small (Mbh — > 0). In principle, by properly tuning the parameter P we can make it as 
small as wanted. Recall that now the solutions have CSS. It should be noted that due 
the the gravitational interaction between the collapsing fluid and the thin shell, we have 
Mbh 7^ m BH + mS BH l i un l ess 6 = 2, which corresponds to null dust. In the latter case, it can 
be shown that by choosing t> = 1 we can make the thin shell disappear, and the collapse is 
purely due to the null fluid. Like the cases with thin shell, by properly tuning the parameter 
P we can make black holes with infinitesimal mass. When £ = 1/2 or 1, which corresponds, 
respectively, to the massless scalar field or to radiation fluid, the solutions reduce to the 



ones considered in 10 



Note that although the mass of black holes takes a scaling form in terms of P, the 
exponent 7 is not uniquely defined. This is because in the present case the solution with P = 
P* = separates black holes from white holes, and the latter is not the result of gravitational 
collapse. Thus, the solutions considered here do not really represent the critical collapse. As 
a result, we can replace P by any function P(P), and for each of such replacements, we will 



have a different 7 ||18|| . However, such replacements do not change the fact that by properly 
tuning the parameter we can make black holes with masses as small as wanted. 



1 While a unique definition of the total mass of a thin shell is still absent, here we simply define it 
as m^fj- 1 = 4irR 2 a. Certainly we can equally use other definitions, such as m s ^fj l = Mbh — m BHi 
but our final conclusions will not depend on them. 
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B. Class B Solutions 



In this case, the first integral of Eq.(pO|) yields 

v = cosh(2r ) - sinh(2r ) tanh(t + t x ), (25) 

where ti is an integration constant. At the moment t = t^H the whole ball collapses inside 
the apparent horizon, and the contribution of the fluid and the shell to the total mass of the 
just formed black holes are given, respectively, by 

m BH = m AH( T AH) = - sinh ?+1 (2r ), 

From the above expressions we can see that for any given 7*0, rn BH and m s BB l are always 
finite and non-zero. Thus, in the present case black holes start to form with a mass gap. It 
should be noted that although m BH is positive, m s BB l is negative. This undesirable feature 
makes the model very unphysical. One may look for other possible junctions. However, 
since the fluid is co-moving, one can always make such junctions across the hypersurface 
r = tq = Const.. Then, the contribution of the collapsing fluid to the total mass of black 
holes will be still given by Eq.(^), and the total mass of the formed balck hole then can be 
written in the form, 

M BH = isinh« +2 (2r ) + M r B f, 

where M r BB denotes the mass contribution of the rest part of the spacetime, which is non- 
negative for any physically reasonable junction. Therefore, in this case for any physical 
junction, the total mass of black holes will be finite and non-zero. 

IV. CONCLUSION 

In this Letter, we have studied two classes of solutions to the Einstein field equations, 
which represent the spherical gravitational collapse of a packet of perfect fluid, accompanied 
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usually by a thin matter shell. The first class of solutions has CSS, and black holes always 
start to form with a zero-mass, while the second class has neither CSS nor DSS, and the 
formation of black holes always starts with a mass gap. The existence of the matter shell 
does not affect our main conclusions. These solutions provide further evidences to support 
our speculation that the formation of black holes always starts with zero-mass for the collapse 
with self- similarities, CSS or DSS. 

It should be noted that none of these two classes of solutions given above represent 
critical collapse. Thus, whether the formation of black holes starts with zero mass or not is 
closely related to the symmetries of the collapse (CSS or DSS), rather than to whether the 
collapse is critical or not. 
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